Bound state formation and nature of the excitonic insulator phase in the extended 

Falicov-Kimball model 
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Motivated by the possibility of pressure-induced exciton condensation in intermediate-valence 
Tm[Se,Te] compounds we study the Falicov-Kimball model extended by a finite f-hole valence band- 
width. Calculating the Frenkel-type exciton propagator we obtain excitonic bound states above 
a characteristic value of the local interband Coulomb attraction. Depending on the system pa- 
rameters coherence between c- and f-states may be established at low temperatures, leading to an 
excitonic insulator phase. We find strong evidence that the excitonic insulator typifies either a BCS 
condensate of electron-hole pairs (weak-coupling regime) or a Bose-Einstein condensate (BEC) of 
preformed excitons (strong-coupling regime), which points towards a BCS-BEC transition scenario 
as Coulomb correlations increase. 

PACS numbers: 71.28.+d, 71.35.-y, 71.35.Lk, 71.30.+h, 71.28.+d. 71.27.+a 
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That excitons in solids might condense into a macro- 
scopic phase-coherent quantum state — the excitonic 
insulator — was theoretically proposed about more than 
four decades ago^, for a recent review see Ref. 0. The 
experimental confirmation has proved challenging, be- 
cause excitonic quasiparticles are not the ground state 
but bound electron-hole excitations that tend to decay 
on a very short timescale. Thus a large number of exci- 
tons has to be created, e.g. by optical pumping, with suf- 
ficiently long lifetimes as a steady-state precondition for 
the Bose-Einstein condensate (BEC) realizing process. 

The obstacles to produce a BEC out of the far- 
off-equilibrium situation caused by optical excitation 
might be circumvented by pressure-induced generation 
of excitons. Hints that pressure-sensitive, narrow-gap 
semiconducting materials, such as intermediate-valent 
TmSeo.45Teo.55, might host an excitonic BEC in solids 
came from a series of electric and thermal transport 
measurements. 3 Fine-tuning the excitonic level, by ap- 
plying pressure, to the level of electrons in the narrow 
4f-valence band, excitons can form near the semiconduc- 
tor semimetal transition in thermodynamical equilibrium 
and might give rise to collective excitonic phases. A 
phase diagram has been deduced out of the resistivity, 
thermal diffusity and heat conductivity data, which con- 
tains, below 20 K and in the pressure range between 5 
and 11 kbar, a superfiuid Bose condensed stated 

The experimental claims for excitonic condensation in 
TmSeo.45Teo.55 have been analysed from a theoretical 
point of view££ Adapting the standard effective-mass, 
(statically) screened Coulomb interaction model to the 
Tm[Se,Te] electron-hole system, the valence-band-hole 
conduction-band-electron mass asymmetry was found 
to suppress the excitonic insulator (EI) phase on the 
semimetallic side, as observed experimentally. But also 
on the semiconducting side, the EI instability might be 
prevented — within this model — by either electron-hole 
liquid phases 6,7 or, at very large electron-hole mass ratios 



(> 100), by Coulomb crystallization.— The effective- mass 
Mott-Wannier-type exciton model neglects, however, im- 
portant band structure effects, intervalley-scattering of 
excitons, as well as exciton-phonon scattering . More- 
over, the excitons in Tm[Se,Te] are rather small-to- 
intermediate sized bound objects (otherwise the exper- 
imentally estimated exciton density of about 1.3 x 10 21 
cm~ 3 would lead to a strongoverlap of the exciton wave 
functions, see Refs. @ and 0). Hence the usual Mott- 
Wannier exciton description seems to be inadequate. 

The onset of an EI phase was invoked quite recently in 
the transition-metal dichalcogenide lT-TiSe2 as driving 
force for the charge-density-wave (CDW) transition^ 

The perhaps minimal lattice model capable of describ- 
ing the generic two-band situation in materials being pos- 
sible candidates for an EI scenario might be the Falicov- 
Kimball model (FKM) , introduced about 40 years ago in 
order to explain the metal-insulator transition in certain 
transition-metal and rare-earth oxidesi 10 ' 11 In its original 
form the model introduces two types of fermions: itiner- 
ant c (or d) electrons and localized f electrons with orbital 
energies e c and ef, respectively. The on-site Coulomb in- 
teraction between c- and f-electrons determines the dis- 
tribution of electrons between these "sub-systems" , and 
therefore may drive a valence transition as observed, e.g., 
in heavy-fermion compounds. To be a good model of the 
mixed- valence state, however, one should build in a co- 
herence between c and f particles J^ This can be achieved 
by a c-f hybridization term V. Alternatively, a finite 
f-bandwidth, which is certainly more realistic than en- 
tirely localized f-electrons, can also induce c-f coherence. 
The FKM with direct f-f hopping is sometimes called ex- 
tended Falicov-Kimball model (EFKM). 

Most notably, it has been suggested that a novel fer- 
roelectric state could be present in the mixed-valence 
phase of the FKM with hybridization.— The origin is 
a non- vanishing (c^/) expectation value, causing a finite 
electrical polarization. In the limit V — > there is no fer- 



roelectric ground state, as was shown in Refs. [l3| and [lj 
in contrast to the findings in Ref. [l2|. Afterwards it has 
been demonstrated that spontaneous electronic ferroelec- 
tricity also exists in the EFKM, provided that the c- and 
f-bands involved have different parity^ 

By means of constrained path Monte Carlo techniques 
the (T = 0) quantum phase diagram of the EFKM 
was calculated in the intermediate-coupling regime for 
one- and two-dimensional (2D) systems, confirming the 
existence of a ferroelectric phase.! 6 A more recent 2D 
Hartree-Fock phase diagram of the EFKMi£ was found 
to agree surprisingly well with the Monte Carlo data, 
supporting mean-field approaches to the 3D EFK M 17 ' 18 
(on the other hand the Hartree-Fock results have been 
questioned by a slave-boson treatment^) . The ferroelec- 
tric state of the EFKM can be viewed as an excitonic 
condensate ((c* f) is an excitonic expectation value since 
/ creates a f-band hole, i.e., the phase with non- vanishing 
polarization is in fact an EI phase). 

Therefore, in this paper, we study the formation of 
excitonic bound states and the nature of the EI phase 
in the framework of the (spinless) 3D EFKM. It can be 
written as 
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where f- and c-orbitals are labeled by the pseudospin vari- 
able a =T,| (or a = ±) with n ka = a\ a a k(T , a fcT = f k , 
and a ki = c k . In Eq. (P), e ka = e a + t a ^ k - /i, 
Ik = \ (cos k x + cos ky + cos k z ) , and /j, is the chemical po- 
tential. The signs of the transfer integrals t a determine 
the type of the gap for large enough \e\— £j.| and/or U and 
of the electronic insulator (ferroelectric) state. Provided 
that tit i < [^t i > 0] we have a direct [indirect] gap 
and the possibility of ferroelectricity (FE) [antiferroelec- 
tricity (AFE)] with ordering vector Q = [Q — (ir, ir, ir)]. 
In what follows, we put e-f = 0, ej > 0, £j_ < 0, if > 0, 
and consider the half-filled band case ^2 n a = 1, where 
n a = (n ia ) = jjY,k( n kv)- 

First we investigate the existence of excitonic bound 
states in the phase without long-range order. To this 
end, we define the creation operator of a Frenkel-type 
exciton by 
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and the exciton commutator Green function 

G x {q,w) = {{b q ;b\)) u . (3) 

In the fermion representation of spins, we have b\ = ST , 
b\ = SZ q = (S+)+, so that Gx(q,w) = ((S+;SZ q )\ u 
may be considered as the negative dynamic pseudospin 
susceptibility. Therefore, to obtain the exciton propaga- 
tor Gx j the calculation of the dynamic spin susceptibility 
in the Hubbard models can be adopted. Taking the ran- 



dom phase approximation we get 

G x (q,uj 
where 
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and to k (q) = s k + q [ — Eki describes the continuum of 
electron-hole excitations. Here, e ka = e ka + Uri- C with 
(nka) = f{e k a) and /(e) = [e e / T + l]~ l . Note that 
G x may be also obtained by the use of the Hamiltonian 
Hi ) = ^2 ka e ka n ka resulting from the Hartree decou- 
pling of the interaction term in Eq. {1} . 

The exciton binding energy is obtained from the poles 
of Gx(q,^) outside the continuum, i.e. from 



,(0) 
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with < lu < u>h{q)\min- Let us consider excitons 
with q = 0. Then, we look for poles with lu = lo x < 

Wfc(0)| m 



= Eg , where the gap E g is given by 
E^ =e i -\t i \-ti + U{ni-n i ). 



(7) 



Here, we are mainly interested in the critical Coulomb 
attraction Ux(T) for the formation of excitonic bound 
states. Eq. ^ is solved numerically for q = 0, using both 
the semielliptic model density of states (DOS) p(e) = 
7^j[l-(£/|*ll) 2 ] 1/2 ( 2|*4.| gives the bandwidth W{) and 
the tight-binding DOS for the 3D cubic lattice (see right 
inset of Fig.[l|). 

In Fig. [T] the boundary Ux(T) for exciton formation 
at U > Ux{T) for both DOS models is plotted (circles, 
diamonds). At low temperatures we find the boundary 
to be rather sensitive to the shape of the DOS, whereas 
at higher temperatures both Ux{T) curves merge. Con- 
sidering a fixed value of U > Ux (T — 0) , with increasing 
temperature the excitons gradually dissociate into sin- 
gle holes and electrons, where at T > Tx(U) the bound 
states are lost. From the cusp of the boundary at the 
point (T S ,U S ) we may suggest an instability against a 
homogeneous phase with long-range order at T < T s . 

The T — Hartree-Fock ground-state phase diagram 
of the 3D half-filled EFKMil^ exhibits— besides full 
f-band and c-band insulator regions at large splittings 
S — \s c — £f\ — two symmetry-broken states: the antici- 
pated excitonic insulator and a CDW. While the CDW 
ground state is stable for all ratios tf/t c at S = 0, it 
becomes rapidly suppressed for 5 > 0, especially if the 
c- and f-bandwidths are comparable.— Since we are in- 
terested in the (homogeneous) condensed excitonic phase 
only, we adjust the parameters S, tf, and U accordingly. 
To model the intermediate-valence situation we choose 
S/\t c \ — 0.4 and tf/\t c \ = 0.8. The almost perfect agree- 
ment between the (2D) Hartree-Fock and path Monte 



Carlo phase diagrams for intermediate couplings 17 might 
justify the application of the Hartree-Fock approach to 
values of the Coulomb attraction U of the order of the 
bandwidth. 

To make contact with previous Hartree-Fock ap- 
proache a 17 ' 18 , we use the equation of motion method for 
the anticommutator Green functions 20 ((a kcr ; a klJ ))ui and 
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the description of the FE EI phase by the order param- 
eter 
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Then, for A ^ 0, we get the gap equation 

N ^ 2E k 

ka 



(9) 
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Figure [T] shows the finite-temperature phase boundary 
of the EI phase obtained by the self-consistent solution of 
the Hartree-Fock Eqs. (I8|)- (fl0|) . In comparison with the 
semielliptic DOS (solid line), the use of the more real- 
istic tight-binding DOS (dashed line) yields a shrinking 
of the EI phase, which corresponds to the behavior of 
the boundary Ux(T) for exciton formation. For U > U s 
we obtain the phase boundary T C (U) coinciding with the 
boundary Tx(U) for exciton formation. This result gives 
a strong argument for the BEC of preformed (tightly 
bound) excitons at T C (U > U s ). On the other hand, 
for U < U s there are no preformed excitons above T c , 
and a BCS-like condensation at T C (U < U s ) takes place, 
i.e., the pair formation and condensation occurs simul- 
taneously. Although the gap equation captures the BCS 
and BEC situation at weak and strong couplings^!, it 
cannot discriminate between them. 

Thus, the existence or non-existence of bound states 
above T c gives strong evidence for a BEC or BCS tran- 
sition scenario at T c , respectively. Moreover, within the 
EI phase, a crossover from a strong-coupling BEC to a 
weak-coupling BCS condensate of electron-hole pairs is 
strongly suggested. 

To describe qualitatively this BEC-BCS crossover re- 
gion, we consider the gap boundary U g (T) (thin dot- 
ted and dashed-dotted lines) resulting from Eg = 
[Eq. 0], where the gap opens for U > U g . Interestingly, 
at T s = T C (U S ) we obtain U g {T s ) = U X {T S ). That is, 
at this point the opening of the gap is accompanied with 
the formation of bound states, whereas for T > T c , U g (T) 
is slightly smaller than Ux(T). From this result we may 
get a crude estimate of the BEC-BCS crossover region by 




FIG. 1: (color online) Phase boundary of the EI phase in the 
EFKM which typifies either a BEC or BCS condensate. Red 
solid (blue dashed) lines are obtained using the semielliptic 
(3D simple cubic) DOS shown in the right inset. The symbols 
give the critical {/-values for exciton formation. Note that 
preformed pairs may exist in the normal phase of the EFKM. 
Black dotted (green dashed-dotted) lines mark the opening of 
the gap Eg (assuming A = 0). For further explanation see 
text. The left inset shows exemplarily the suppression of the 
order parameter A with increasing temperature on the track 
from C — > A, i.e. at fixed U, where a second-order phase 
transition is obtained. 



extrapolating U g (T) into the EI phase. Solving Eq. ^ 
at a fixed T < T c , for U in the region U g < U < U x 
we get negative pole energies lox which indicates the in- 
stability of the normal phase with bound states against 
the long-range ordered EI phase. Moreover, for U <U g 
no solution can be found which may be indicative for 
an instability towards a BCS-type EI state. Thus, the 
BEC-BCS crossover in the EI phase should occur in the 
neighborhood of the U g (T) line. 

In comparison to the phase boundary obtained within 
the simple effective-mass, Mott-Wannier-type models 
the EI phase of the EFKM is confined at zero temper- 
ature on the weak-coupling side, because of the finite 
f- and c-bandwidths. While the shape of the EI dome 
approximates the Tm[Te,Se] phase diagram constructed 
from the experimental data, the absolute transition tem- 
peratures are overestimated, of course, by any mean-field 
approach. The homogeneous EI phase shrinks as 5 as tf 
becomes smaller at fixed S, but it does not disappear M- 

Figure gives the partial f- and c-electron DOS at 
various characteristic points A-F of the phase diagram 
shown in Fig. [1] The high-temperature phase may be 
viewed as a metal/semimetal (panel A) or a semiconduc- 
tor (panel B) in the weak- or intermediate-to-strong in- 
teraction regime, respectively. The EI phase shows com- 
pletely different behavior. As can be seen from panel C, 
a correlation-induced "hybridization" gap opens, indicat- 
ing long-range order (non- vanishing f-c-polarization) . As 
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FIG. 2: (color online) DOS for f-band ([black] solid curves) 
and c-band electrons ([red] dashed curves) at the points 
marked by A to F in the T-U plane in Fig.[T] Band structure 
parameters are 8 — 0.4|i c |, tf = 0.8|£ c |. For comparison, the 
results for 5 — 0.1|t c | an d tf = 0.05|t c | ([green] double-dot- 
dashed and [blue] dot-dashed curves) have been included in 
panels E and F; here a CDW will become the Hartree-Fock 
ground state as U increases from E to F. All data are obtained 
for the semielliptic DOS and total filling n = Uf +n c = 1; the 
arrows mark the positions of the chemical potential. 



the temperature increases the gap weakens and finally 
closes at T — T c . The pronounced c-f-state mixing and 
strong enhancement of the DOS at the upper/lower va- 
lence/conducting band edges is reminiscent of a BCS-like 
structure evolving from a (semi-) metallic state with a 



large Fermi surface above T c (see panel A) . This may be 
in favor of a BCS pairing in the weak-coupling region of 
the EI phase, as discussed above. By contrast the DOS 
shown in panel D clearly evolves from a gapped high- 
temperature phase. Finally, in panel E [F] the partial f- 
and c-electron DOS at T ~ below [above] the EI phase 
are depicted, where the system behaves as a metal or 
semimetal [band insulator or semiconductor]. Note that 
the splitting of c- and f-bands in panel F is not caused 
by 6 (being the same as in E), but is due to the Hartree 
shift oc Uijif — n c ). 

To summarize, in this work, we attempted to link ex- 
perimental hints for excitonic condensation to recent the- 
oretical studies of electronic ferroelectricity in the ex- 
tended Falicov-Kimball model. We analyzed the finite- 
temperature phase diagram and argued that a finite 
f-bandwidth in combination with a short-range inter- 
band Coulomb attraction between (heavy) valence-band 
holes and (light) conduction-band electrons may lead 
to f-c-band coherence and an excitonic insulator low- 
temperature phase. Most noteworthy we established the 
existence of excitonic bound states for the EFKM on the 
semiconductor side of the semimetal semiconductor tran- 
sition above T c , and suggested a BCS-BEC crossover sce- 
nario within the condensed state. As a consequence, we 
expect pronounced transport anomalies in the transition 
regime both in the low- and high-temperature phases, 
which should be studied in the framework of the EFKM 
in future work. 
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